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' . 1 Introduction 

Oh 

D . __ 

This talk is intended as a small advert for using holography [ 1 2 1 to study strongly coupled systems 
such as QCD. Holography provides new insights into how to model non-abelian gauge theory in en- 
vironments that were previously beyond our reach. To stress that power I will discuss some examples 
from my own work [3 4 5 6 1 . Firstly I will show that holography allows computation at finite chemical 
potential and compute the phase diagram of a range of strongly coupled systems. Secondly I shall look 
at time dependent problems during finite temperature phase transitions. The rigorous examples I shall 
present are not QCD but other strongly coupled gauge theories (although still with gauge fields and 
quarks). The ideas behind holography apply equally well to QCD although we do not have a rigorous 
derivation of the dual - one can nevertheless model build sensibly using the formalism. 

2 Gauge/Gravity Duality 



Holography [ 1 2| provides a unique way to model strongly coupled gauge theories. The renormaliza- 
tion group scale of the gauge theory is treated as a space-time dimension. The conformal symmetry of 
classical gauge theory is realized as a symmetry of the AdSs metric 



R 

ds 2 = ^dx 2 M + —dr 2 , (1) 



where R is the radius of curvature of the space. We can think of the space as a box with r corresponding 
to energy scale. At any fixed r we see a 3 + Id theory living in the X3 + i directions - large r is the UV 
of theory whilst small r is the IR. The dilatation symmetries in the classical gauge theory act on the 
space and the fields 

x -> e a x, A* -> e- a A» (2) 
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Fig. 1. Quarks can be added to N = 4 SYM as 3-7 strings in the D3-D7 system. 



but are realized on the AdS space as 



r — > e 



r. 



(3) 



Fields in the AdS space represent gauge invariant operators and sources in the gauge theory such as 
the gauge coupling g 2 YU or TrF 2 . 

The most clear cut example of this formalism is for the conformal N = 4 supersymmetric gauge 
theory with a large number of colours (the theory has a gauge field, 6 real scalars and 4 two-component 
gauginos all in the adjoint representation of SU(N)). In string theory this theory lives on the world- 
volume of a stack of D3 branes in the limit where the string tension is taken to infinity. The same D3 
branes' tensions source an AdSs space in the ten dimensional string theory. When there are a large 
number of D3s (N — > oo) the radius of the geometry R — > oo relative to the string length and a classical 
gravitational description emerges. Maldacena fl] made the remarkable leap of faith that the surface 
gauge theory and linked closed string theory in AdSs were alternative descriptions of each other. The 
full power of this statement is that the weakly coupled degrees of freedom in the AdS space (strictly 
supergravity modes in the tension goes to infinity limit) and their interactions are precisely known 
from the string theory and describe the strongly coupled gauge dynamics. 

Since the original insight a great deal of work has been done in string theory. It is now believed 
that the holographic principle works for a wide range of gauge theories [7] including all so called 
deformations |8 9 10] of N — 4 SYM, but also for any theory realizable in the world-volume of branes. 
In fact it is clear that my first paragraph applies as a starting point for any gauge theory although there 
is no assurance that the neccessary deformations of AdS will leave a weakly curved space, nor that one 
can identify the bulk fields and their interactions. However, holography does provide a new framework 
for constructing effective descriptions of strong dynamics and there is ample evidence that sensible 
quantitative results can be found in this way for QCD 111 11121 . Perhaps the most powerful insight is 
that holography allows computation in environments where other techniques fail. Here we compute at 
finite density and look at time dependent problems to show the technique's versatility. 



Quark fields must be an integral part of any description of QCD. At a phenomenological level we can 
introduce a scalar field, i^, into the AdS space that is dual to the quark mass and {qq) condensate. In 
the UV we can propose the action 



3 Quarks 




(4) 



which for the AdS background leads to solutions 



P(4-7) 




the integration constants have the correct dimensions to play the roles of the mass and condensate. The 
IR though should be characterized by more subtle interactions of the scalar with itself and in principle 
the background geometry. 

Top down string constructions II13I14I15I16I17I181I provide some insight into the form the action 
could take. Quark N — 2 hypermultiplets can be added to the N = 4 SYM theory through the 
introduction of D7 branes as shown in Fig Q] The minimum length of the 37 string is a measure of 
the quark mass (its minimum energy is given by the length times tension). In the quenched Nf <K N c 
limit the D7 branes can be treated as probes in the AdS geometry. The embedding is described by a 
single function L(p) representing the separation of the D3 and D7 in the 89 directions, L, as one moves 
along the D7 in the 4567 directions, p. The action for this embedding is given by the Dirac-Born-Infeld 
action of the brane (in Einstein frame) 

S D1 = ~T 7 f V-detP[GU = j d*x dp pV ^1 + {d p Lf (6) 

where Tj = l/(27r) 7 a 4 and Tj = 2n 2 T-ilg 2 D when we have integrated over the 3-sphere on the D7. 
is the dilaton that is dual to the gauge coupling. Expanding at large p one can identify if/ = L/p. In 
an AdS background the regular solution (satisfying d p L(Q) = 0) is the flat embedding shown in FigQ] 
Massless quarks occur when the D7 lies to intersect the D3 at p — 0. 



4 Chiral Symmetry Breaking and Confinement 

An example of a complete top down D3-D7 model that displays chiral symmetry breaking is provided 
by the imposition of a background magnetic field [19 20 3 |. The D7 DBI action contains a world- 
volume gauge field that is dual to the operators qy^q and its source, a background U(l) baryon number 
gauge field. It is straightforward therefore to study the model in a background B field. In fact the DBI 
action for the D7 that results is that in (0 with an effective un-back-reacted dilaton 

The presence of such a dilaton disfavours D7 embeddings that penetrate to r — y/p 2 + L 2 = since 
the action grows as r — > 0. The D7 brane therefore bends away from the D3 as shown in Fig [2] We 
can naively see that a configuration that at large p (the UV) describes a massless quark is modified at 
small p (the IR) and there in fact is no zero length 37 string. The quark has acquired an effective IR 
mass due to the chiral symmetry breaking (the asymptotic curvature of the embedding corresponds to 
the presence of a quark condensate in the spirit of the second solution in (0 ). 



Chiral symmetry breaking is induced in QCD by the running gauge coupling which breaks con- 
formal invariance, introduces strong coupling and the scale Aqcd- In a phenomenological spirit we 
can introduce a running coupling by imposing a non-trivial profile for the dilaton (which we will not 
backreact on the geometry) - for example we can create a step in the dilaton profile [4] 

e* = A + 1 - A tanh [r(r - A)] (8) 

where A generates the scale equivalent to Aqcd, A is the step size, and r controls the energy regime 
over which the change occurs. The D7 brane dynamics in this background is similar to the background 
B field case - the dilaton grows into the IR and repels the D7 brane again triggering chiral symmetry 
breaking as shown schematically in Fig [2] 

The dilaton also enters into the Nambu Goto action for a string. A quark anti-quark pair and their 
interactions are represented by a 77 string as shown in Fig|2l 1 2112211 . Maldacena |22| has shown that 
in AdS the string droops to deeper values of r the larger the quark anti-quark separation, d, is and the 
string energy scales as l/d consistent with a Coulomb law in the conformal theory. When a dilaton 
profile that causes repulsion of the string from the interior of AdS is introduced, the string tends to lie 
outside the repulsion radius (often called a "wall")- The energy of the string then grows as d as the 
quarks are separated since the length of string on the wall is simply lengthened at constant energy cost 
per unit length. Our phenomenological dilaton profile will therefore lead to confining behaviour too 
0. 



5 Temperature and Density 

It has long been known how to deform the AdS space to describe the N — 4 gauge theory at finite tem- 
perature (2). The correct geometry is an AdS Schwarzschild black hole. The Hawking temperature of 
the black hole effectively heats up the dual and represents the gauge theory temperature. The geometry 
is given by 

dg2 = ^2 {~f dtl + dx l) + Jri dr2 i f = 1 " ^ (9) 

When D7 branes are introduced into this geometry a new possible configuration exists in which the 
D7 brane falls into the black hole horizon [ 17 23 24]. On the gauge theory side this phase represents 
one where the quark degrees of freedom have become unbound and have merged with the background 
gluonic thermal bath 1251261 . 

Baryon number chemical potential and density are also easily introduced using the D7 world- 
volume gauge field previously used to introduce magnetic field [27 28 29 30]. One allows a non-zero 
profile for A, - for example if it is a constant, fi, then the quark kinetic term replicates a chemical 
potential term 

qy°A t q -» qy°nq (10) 

ie a shift in the quark vacuum energy of fi, corresponding to a Fermi surface. 
The relevant gauge field enters the DBI action as 

Sdi = -r 7 J </V ^-det(P[G] ab + F ah ) (11) 

and in the full solutions has UV asymptotic behaviour A, = fi + A . Here d represents the quark 

density operator qy°q. In fact since the action only depends on derivatives of A, there is a conserved 
quantity that one may choose, which is precisely d. If one sets d then one can solve a single ODE for 
the embedding L(p) and then solve for yu(p) using the L(p) solution. This is easier than setting /i and 
solving two coupled ODEs for L(p), A,(p). When the density d is non-zero, the D7 brane embedding 
L(p) is deformed so that it always ends at r = 113011 . ie on the D3 branes, or at the black hole horizon 
if there is also temperature present. The reason is that a non-zero quark density corresponds to strings 
stretching between the D3 and the D7 - in the infinite tension limit they force the D3 and D7 together. 




6 Phase Diagrams 

We now have all the ingredients we need to determine the phase structure of models of gauge theories 
with chiral symmetry breaking [345]. For a given choice of dilaton one can fix T and fi and then 
numerically seek regular embeddings of the D7 brane. We will work with massless quarks so that 
L — > as p — > oo. The configurations that contribute to the phase structure are: the fiat embedding, 
L = 0, which preserves chiral symmetry; configurations that bend off axis ending off the black hole 
horizon so that d p L(Q) = 0, which are chiral symmetry breaking configurations; between these two 
configurations are solutions that curve off axis but end on the blackhole horizon representing a phase 
with chiral symmetry breaking and non-zero quark density. Finding all such embeddings across the 
phase plane is numerical intensive but straightforward. The true vacuum can be found by comparing 
the configurations' free energy, given by minus the DBI action. 

In Fig[3]we display the phase diagram for the N — 4 SYM theory with N — 2 quark hypermulti- 
plets in the presence of a magnetic field (ie we use the dilaton form in GJ|)[3|. It shows all three phases 
discussed separated by a complex mixture of first order (blue lines) and second order (red lines) tran- 
sitions, including two critical points. It is worth stressing that this is an example of what a powerful 
technique holography is - no other method would allow the computation of this phase diagram! 

Note that in this framework to move from a chiral symmetry breaking, d — 0, configuration to the 
chirally symmetric phase smoothly requires an intermediate phase where the D7 curves off axis and 
ends on the black hole. That suggests the presence of a phase with a non-zero quark condensate and 
quark density is very likely in a range of gauge theories. 

I would argue that the B field catalysis of chiral symmetry breaking makes a reasonable first stab 
at QCD-like dynamics. At these transitions QCD is strongly coupled and characterized by the chiral 
symmetry breaking scale Aqcd - here we have a strongly coupled gauge theory with a conformal 
symmetry breaking scale yB. Of course our theory has additional degrees of freedom associated with 
the superpartners but at non-zero T and fi supersymmetry is broken so there are no special cancellations 
present in the theory. The phases and their positions in the plane do indeed match expectations in QCD. 
The details though are not correct - in QCD one expects a first order transition at low T and high fi but 
a second order transition at low fi and high T. 

It is natural to ask how robust the phase structures of the model are. To test this we will move to 
the phenomenological ansatz for the coupling in ^ which allows us to vary the strength of conformal 
symmetry breaking by hand [4]. Of course we are now moving away from a precise description to 
modelling but with many of the key ingredients of the holographic chiral symmetry breaking mech- 



Fig. 4. Plots of three possible phase diagrams for the choices A = 3, 15, 30 |4 1 . Large (small) A gives second (first) 
order transition at low T. f = 1, A = 1.715. 




Fig. 5. Sample phase diagrams for theories with none zero a (a = 2.2, a = 3) [ 4 1 . 



anism intact. The computational methods remain the same In Fig [4] we show three examples of 
the phase structure at fixed A, r but for different values of the step height in the dilaton ansatz (0. 
The three phases correspond in the obvious fashion to those shown in the B field case in Fig [3] The 
main conclusion is that reducing the strength of the conformal symmetry breaking (ie reducing A) 
in the gauge coupling leads to a preference for a first order chiral phase transition. This model cer- 
tainly shows that it is relatively easy to engineer a first order transition with chemical potential at zero 
temperature. 

In fact with a little more phenomenological tinkering a second order transition with temperature 
at low density is also achievable. Away from the probe approximation one should backreact the D7 
branes on the AdS geometry. Some serious attempts to include this backreaction for the N — 2 SYM 
model can be found in [ 14 32 33 34]. Here though in our phenomenological model we simply note that 
the D7 break the symmetry between the p and L directions of the geometry. This symmetry breaking 
is not present in our black hole ansatz in (O. A simple way to introduce it is to write J4) 

r 2 =p 2 + L 2 /a (12) 

in the emblackening factor / of the black hole geometry. By dialling the phenomenological parameter 
a one can break the p, L symmetry. If we now repeat our analysis we find the phase diagrams shown 
in Fig [5] The second order transition at low density in the figure on the right is the main result. This 
sensibly motivated phenomenological tinkering suggests that the phase structure of strongly coupled 
gauge theories is quite sensitive to the precise conformal symmetry breaking structure, but equally that 
it is possible to model a range of phase diagrams using holography. 

A final twist to the story was provided in where configurations corresponding to a non-zero 
density of baryons were included. Witten described baryons by a D5 brane wrapped on the S 5 of the 
dual geometry [35 36] . Such a configuration must have a surface gauge field sourced by N fundamental 
strings ending on it to be consistent, suggesting it can be thought of as a baryon. In [5] new configura- 
tions of a D7 ending on such a wrapped D5 (the fundamental strings pull the D5 and D7 to meet at a 
point, subject to a force balancing condition |37 38 39]) were included in the model so far discussed. 
For particular choices of parameter space a QCD-like baryonic phase was achieved as shown in Fig|6] 
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Fig. 6. Sample phase diagram including a baryonic density phase 0. 



The transition to the phase with baryon density is second order whilst in QCD it is known to be first 
order at low T - the difference is due to the absence of inter-baryonic interactions in our description 
(attempts to include such interactions can be found in [40 1). 

In conclusion to this section then we see that holography allows the full computation of the phase 
diagrams of strongly coupled gauge theories. The case of the N — 2 theory with a magnetic field is 
believed to be a complete computation. Beyond that theory, we have studied how sensitive the phase 
structure is to phenomenological variation of the conformal symmetry breaking. Broadly the phases 
present are very robust (and match those expected in QCD). The orders of the transitions are much 
more sensitive but this does then allow us to cook models of the QCD phase diagram. Further we learn 
that the weaker the conformal symmetry breaking of a theory the more likely the chiral restoration 
transition is to be first order. 



7 Out of Equilibrium Dynamics 

Finally I would like to discuss an additional boon of holographic modelling. It is now relatively 
easy, using holography, to study time dependent problems in strongly coupled theories since they 
are mapped to time dependence in weakly coupled classical systems. As a simple example of the types 
of computation that are now open to study in [ 6 ] we looked at the thermal, first order, chiral restoration 
phase transition in the N — 2 gauge theory with a magnetic field. 

We were able to study this system's time evolution thanks to the work in 114 1142143 II where a time 
dependent geometry describing a growing/shrinking black hole in AdS was computed. This geometry 
has a black hole whose horizon moves away from the boundary in time as the N — 4 plasma it 
describes expands and cools. The time reversed solution describes a heating plasma and I will find it 
useful to discuss that scenario too below. The geometry is a late time expansion (when the black hole 
is small) in powers of inverse time. By controlling the strength of the B field on the D7 probes we can 
arrange to place the chiral phase transition at any point in the evolution so the expansion is sufficient 
to fully study the transition. 

The geometry is given by 




(13) 



where z— 1/r. At late times, the coefficients can be expanded to first order as 




(14) 



c(t,z) = \n(l + v 4 /3) + (-2r ]0 £ 




Fig. 7. Results for time dependence of the condensate through the first order phase transition in the N = 2 gauge 
theory with magnetic field. Left: The condensate c{t) in the equilibrium and non-equilibrium (blue) approaches. 
In the equilibrium description, the different branches correspond to the extrema of the potential V(c) on the Right. 
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with 



_ z 1/4 
V =Jj3 e ' 
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^ - 2 l/2 3 3/4 ■ (15 ^ 

€q is a free parameter of mass dimension 8/3 and is related to the energy density, while 770 is related 
to the shear viscosity, v is a scaling parameter valid at large t. Note that a(t,z),b(t,z) and c(t, z) are 
expanded around t — 00 in powers of 1/f 2 ^ 3 as 
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(16) 



Expanding g„ near the boundary (z = 0) we can recover a time-dependent emblackening factor, 
which describes a moving horizon. The size of the horizon determines the time-dependence of the 
'temperature' as 1421441 
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The DBI action for the D7 branes schematically reads [6 1 



J d s ^-det(P[G] ab ) ~ J dtdpp^A^jl +(d p L) 2 



(p 2 + L 2 ) 2 



(18) 



with A, B complicated but computable functions of p, t. The evolution of the D7 embedding is now 
given by a 1 + 1 dimensional PDE which can be solved by standard numerical packages. 

The simplest example of a quantity one can study is the quark condensate (c) as a function of time 
as the plasma heats or cools. A first approximation is simply to take the expression for T(t) from ( [T71 > 
and for each time use the equilibrium computation for the condensate (ie use a static black hole of that 
temperature, find the D7 embedding and read off the condensate from the asymptotic curvature). That 
computation is represented by the black curve in the left figure of Fig|7] The time axis can be read from 
left to right for a cooling plasma or right to left for a heating plasma. At large time the preferred vacuum 
has a non-zero condensate. Around t = 25 the first order transition occurs to the chirally symmetric 
condensate. The behaviour is consistent with the schematic temperature dependence of the potential 
shown in the right hand figure in Fig [7] At large times c — is a maximum of the potential which 
has a minimum at non-zero c. Around the transition there are two minima one with c = and one 
with it non-zero, and a maximum between. At low temperatures the non-zero c minimum disappears 
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Fig. 8. Plots of the IR position of the D7 brane against time for a number of large time initial conditions with 
different energy (6) . The thick black line is the horizon. 

and c — is the unique vacuum. This is a standard first order phase transition with c jumping from a 
non-zero value to zero around the transition point. 

Of course in a real system c does not jump discontinuously. The transition proceeds by bubble 
formation. The full PDE solutions can be used to study this process. For example if one begins at very 
large time with the static vacuum configuration at T — one can watch the evolution as the vacuum 
is heated. The blue curve in the left hand figure of Fig Q shows the result. It follows the equilibrium 
expectations until the vacuum at non-zero c ceases to exist at which point it smoothly transitions to 
the c — configuration. This represents a super cooled transition. There is no thermal excitation of the 
vacuum about the minimum at large time, the vacuum sits in the right hand minimum of the potential 
shown on the right of Fig|7]until the temperature where that minimum ceases to exist when it rolls to 
c = 0. 

One can add in some thermal excitation by giving the D7 initial condition a small amount of 
non-zero L - at large times it then oscillates about the true vacuum embedding (performs harmonic 
oscillations about the minimum of the potential at non-zero c in Fig |7). Depending upon how much 
thermal energy it has, the configuration will move over the hump in the potential near the transition 
point earlier or later. This is a simplistic description of bubble formation in hotter regions. In Fig[8]we 
plot L(0) for D7 configurations with varying initial thermal energy and we indeed see that they hit the 
black hole horizon earlier the more energy they have ie they escape to the c — vacuum earlier in their 
evolution. 

In principle one could study systems with spatially dependent initial conditions for the embedding 
and one would then observe full bubble formation, bubble collapse before the transition point, and 
finally bubble growth to drive the transition. These are the phenomena that are familiar from watching 
a kettle boil of course, but here, remarkably, we are computing them for a strongly coupled gauge 
plasma. 



8 Conclusions 

Holography represents a new formalism in which to present renormalization group flow in strongly 
coupled gauge theories and, at least for a class of gauge theories derivable from string theory, is 
known to be a weakly coupled description. It is straightforward to compute in time dependent problems 
and in environments with finite chemical potential that are a major obstacle for other approaches. As 
examples of the technique's power, I have presented models of the phase diagrams for QCD, and more 
rigorous computations of an N — 2 gauge theory in the presence of a magnetic field, including studies 
of the first order thermal transition which proceeds via bubble formation. 
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